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Abstract 

This  paper  describes  a  computer  modeling  and  simulation  system 
that  supports  computational  steering,  which  is  an  effort  to  make 
the  typical  simulation  workflow  more  efficient.  Our  system  pro¬ 
vides  an  interface  that  allows  scientists  to  perform  all  of  the  steps 
in  the  simulation  process  in  parallel  and  online.  It  uses  a  standard 
network  flow  visualization  package,  which  has  been  extended  to 
display  graphical  output  in  an  immersive  virtual  environment  such 
as  a  CAVE.  Our  system  allows  scientists  to  interactively  manipu¬ 
late  simulation  parameters  and  observe  the  results.  It  also  supports 
inverse  steering,  where  the  user  specifies  the  desired  simulation 
result,  and  the  system  searches  for  the  simulation  parameters  that 
achieve  this  result.  Taken  together,  these  capabilities  allow  scien¬ 
tists  to  more  efficiently  and  effectively  understand  model  behavior, 
as  well  as  to  search  through  simulation  parameter  space. 

This  paper  is  also  a  case  study  of  applying  our  system  to  the 
problem  of  simulating  microwave  interactions  with  missile  bodies. 
Because  these  interactions  are  difficult  to  study  experimentally,  and 
have  important  effects  on  missile  electronics,  there  is  a  strong  desire 
to  develop  and  validate  simulation  models  of  this  phenomena. 
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1  Introduction 

This  paper  describes  computational  steering,  which  is  an  effort  to 
make  the  typical  modeling  and  simulation  workflow  more  efficient. 
The  particular  modeling  and  simulation  application  where  we  have 
applied  our  efforts  is  the  interaction  of  high-power  microwaves  with 
missile  bodies.  This  application  has  the  goals  of  1)  better  under¬ 
standing  this  complicated  and  hard-to-measure  phenomenon,  and  2) 
driving  the  design  of  both  better  missiles  and  potential  anti-missile 
technologies. 

1.1  Computational  Steering 

Computational  steering  is  an  evolution  of  the  traditional  batch  se¬ 
quence  —  1)  execute  the  model,  2)  analyze  the  results,  3)  mod¬ 
ify  the  parameters  —  that  computational  scientists  have  engaged  in 
for  decades.  This  traditional  batch  process  is  not  always  effective, 
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especially  when  executing  the  model  is  very  time-consuming.  In 
computational  steering  the  same  three  steps  all  take  place  more  or 
less  concurrently  and  online.  This  places  the  scientist  in  the  middle 
of  a  running  experiment,  and  allows  them  to  steer  the  parameters 
while  the  experiment  is  underway.  The  ultimate  goal  is  to  enable 
the  scientist  to  achieve  the  objectives  of  the  experiment  more  effi¬ 
ciently  and  effectively.  An  overview  of  computational  steering  is 
given  by  Parker  et  al.  [7], 

Computational  steering  can  be  done  in  two  directions.  In  for¬ 
ward  steering,  the  scientist  directly  manipulates  the  model  steering 
parameters,  based  on  the  observed  behavior  of  the  model.  If  the 
model  execution  time  is  sufficiently  small,  and  the  graphical  visu¬ 
alization  sufficiently  revealing,  this  can  allow  the  scientist  to  effi¬ 
ciently  explore  the  steering  parameter  space,  and  discover  parame¬ 
ter  settings  that  are  optimal  in  some  sense.  However,  when  model 
behavior  is  very  complex,  when  there  are  many  steering  parameters, 
or  when  it  takes  too  long  to  execute  the  model,  it  becomes  infeasible 
for  the  scientist  to  manually  steer  to  the  desired  solution.  For  these 
situations  inverse  steering  is  needed.  With  inverse  steering,  the  sci¬ 
entist  defines  desired  model  behavior,  and  the  system  automatically 
searches  the  steering  parameter  space,  looking  for  parameter  set¬ 
tings  that  produce  the  desired  behavior.  This  effectively  exchanges 
the  input  and  output  of  the  forward  steering  process. 

1.2  Application  Area 

In  this  work,  we  are  simulating  the  interaction  of  high-power  mi¬ 
crowaves  with  missile  body  casings.  These  interactions  have  po¬ 
tentially  important  consequences  for  the  electronics  of  the  missile. 
High  power  microwave  radiation  could  come  from  either  a  nearby 
nuclear  explosion  or  an  anti-missile  weapon  system.  Microwaves 
are  just  electromagnetic  (EM)  radiation  in  a  certain  frequency  band. 
In  our  simulations,  we  set  up  a  simplified  model  of  a  missile.  We 
then  define  an  EM  radiation  source,  and  begin  the  simulation  such 
that  at  time  step  0,  radiation  is  just  beginning  to  arrive  at  the  model. 
The  radiation  encounters  the  model  and  begins  a  complex  series  of 
reflections.  Typically,  after  several  hundred  time  steps,  the  simula¬ 
tion  approaches  a  steady  state,  where,  in  the  vicinity  of  the  model, 
the  EM  field  does  not  change  with  additional  time  steps,  except  for 
an  oscillation  inherent  in  the  frequency  of  the  EM  radiation  source. 
The  simulation  typically  reveals  hot  zones,  which  are  regions  of 
local  maximum  EM  energy.  The  locations  and  intensities  of  these 
hot  zones  are  an  important  application  design  parameter  which  our 
work  aims  to  supply.  With  this  knowledge,  engineers  can  design  the 
missile  so  that  critical  electronics  are  protected  from  the  hot  zones. 

2  System  Design 

2.1  System  Architecture 

Figure  1  illustrates  the  basic  architecture  of  our  system.  The  sci¬ 
entist  interacts  with  the  system  using  a  desktop  system  or  a  CAVE. 
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Figure  1:  Overall  system  architecture. 


The  steering  control  and  display  manager  contains  the  interface 
between  the  scientist  and  the  rest  of  the  system  components.  This 
process  is  implemented  within  AVS,  a  standard  pipeline-based  sci¬ 
entific  visualization  package.  Color  plate  1  shows  an  image  of 
the  AVS  control  panel.  Model  execution  interfaces  with  the  other 
system  components,  providing  both  user  control  and  observation. 
The  scientist  controls  the  model  execution  —  starting,  stopping, 
parameter  adjustment,  etc. —  directly  from  the  AVS  control  panel. 
The  parameter  steering  interface  converts  these  control  parame¬ 
ter  adjustments  into  appropriate  model  updates  within  the  model 
database.  These  updates  could  be  as  simple  as  setting  a  single 
variable,  or  as  complex  as  setting  a  number  of  arrays.  The  objec¬ 
tive  function  seeker  provides  inverse  steering.  It  obtains  the  current 
model  state  from  the  executing  process,  evaluates  that  state  accord¬ 
ing  to  the  desired  objective  function,  and  then  adjusts  control  pa¬ 
rameters  accordingly. 

Our  system  runs  on  a  high-end  Silicon  Graphics  Onyx  2,  with  6 
processors  and  2  graphics  pipes.  This  machine  also  runs  the  CAVE, 
and  we  can  use  it  to  execute  smaller  models  locally.  As  shown  in 
Figure  1,  we  use  a  remote  supercomputer  (described  in  Section  2.2 
below),  for  moderate  to  larger  models.  A  high-speed  ATM  network 
connects  the  two  machines. 

2.2  Computational  Issues 

Simulating  EM  field  interactions  requires  solving  Maxwell’s  equa¬ 
tions  for  the  propagation  of  electromagnetic  fields.  This  has  been 
a  common  simulation  application  for  30  years,  and  many  standard 
algorithms  for  solving  Maxwell’s  equations  exist  [9],  However,  the 
limited  accuracy  of  standard  algorithms  demand  a  dense  computa¬ 
tional  grid,  which  consequently  makes  the  execution  time  unten¬ 
able  for  even  moderately- sized  models.  In  the  early  1990’s  James 
Cole  developed  the  HASP  (High-Accuracy  Scattering  and  Propa¬ 
gation)  algorithm  [1,  2,  3],  at  the  Naval  Research  Lab  (NRL).  The 
HASP  algorithm  is  up  to  10,000  times  more  accurate  than  stan¬ 
dard  approaches  [3].  This  increase  in  accuracy  permits  the  use  of 
a  relatively  coarse  grid,  which  greatly  improves  memory  use  and 
execution  time.  In  addition,  greater  algorithmic  stability  allows  a 
reduction  in  the  number  of  necessary  iterations.  Furthermore,  the 
structure  of  the  algorithm  makes  it  very  amenable  to  parallel  pro¬ 
cessing. 

As  described  in  Section  1.2  above,  we  are  studying  the  effects  of 
high-power  microwaves  on  missile  bodies.  Because  experimental 
studies  of  EM  interactions  with  actual  missile  structures  are  very 
difficult  to  arrange,  cooperating  researchers  have  developed  a  stan¬ 
dard  test  body,  termed  the  standard  cylinder  [6],  The  standard 
cylinder  allows  multiple  research  groups  to  experiment  with  a  com¬ 
mon,  unclassified,  test  body,  and  to  easily  compare  results.  Figure  2 
is  a  schematic  drawing  of  the  standard  cylinder.  When  tested  exper¬ 
imentally,  the  field  strength  probe  measures  the  EM  field  at  various 
places  along  the  center  axis.  This  probe  must  be  included  in  the 
model,  as  it  effects  the  EM  field  inside  the  cylinder.  While  we  have 
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Figure  2:  The  standard  cylinder,  a  commonly-used  model  for  test¬ 
ing  microwave  interactions  with  missile  bodies. 


used  a  variety  of  models  in  our  system,  we  have  most  extensively 
studied  the  standard  cylinder  model.  All  of  the  Color  Plate  images 
in  this  paper  are  of  the  standard  cylinder. 

The  efficiency  of  the  HASP  algorithm  allows  the  work  reported 
here.  To  achieve  acceptable  computation  times,  we  implemented 
the  HASP  algorithm  on  NRL’s  128  processor  Silicon  Graphics  Ori¬ 
gin  2000  supercomputer,  using  both  the  automatic  parallelization 
option  of  the  SGI  MISPRO  Fortran90  compiler,  as  well  as  aggres¬ 
sive  optimization  techniques  to  speed  the  scalar  code  running  on 
each  processor.  For  the  standard  cylinder,  this  implementation  cal¬ 
culates  each  time  step  in  a  few  tenths  of  a  second,  and  a  steady  state 
result  (~  10  billion  operations)  in  a  few  minutes. 

2.3  Visualization  Interface 

Color  Plate  1  shows  the  user  interface  for  the  system.  The  control 
panel  to  the  left  controls  both  the  simulation  execution  and  the  steer¬ 
ing  parameters.  Scientists  can  start,  stop,  and  resume  the  simulation 
at  any  time  step,  and  can  save  the  system  state  to  disk.  They  also  use 
this  control  panel  to  modify  the  five  scalar  steering  parameters:  mi¬ 
crowave  source  frequency,  source  direction  angle  polariza¬ 

tion  angle,  and  probe  insertion  distance  (note,  in  Color  Plate  1,  the 
five  dial  controls  in  the  steering  control  panel).  The  control  panel 
on  the  top  controls  various  display  parameters,  such  as  slice  plane 
values  and  glyphs.  The  bottom  of  Color  Plate  1  shows  the  AVS  net¬ 
work  that  controls  receiving  and  viewing  the  results.  To  the  right  are 
two  views  generated  by  the  two  geometry  viewer  modules;  a  pair  of 
blue  arrows  indicate  which  viewer  produces  each  view.  Also  shown 
is  another  module,  termed  the  VR  Viewer ,  which  generates  output 
for  a  CAVE  display.  This  module,  developed  at  NRL  [5],  replaces 
the  AVS  geometry  viewer  and  displays  information  on  any  VR  dis¬ 
play.  The  VR  Viewer  uses  the  OpenGL  Volumizer  library  to  render 
volumes  using  3D  texture  memory,  and  allows  the  scientist  to  inter¬ 
act  with  their  dataset  via  a  tracked,  hand-held  3D  wand.  Color  Plate 
3  shows  a  scientist  interacting  with  the  data  in  our  CAVE  facility.  It 
also  shows  how  we  sometimes  use  one  or  two  CAVE  walls  to  dis¬ 
play  steering  parameters,  color  maps,  and  other  AVS  widgets.  We 
mounted  a  keyboard  and  mouse  on  a  portable  stand  to  allow  users 
to  manipulate  these  interface  elements. 

2.4  Inverse  Steering 

In  general,  the  scientists  and  designers  using  our  system  are  inter¬ 
ested  in  two  different  types  of  phenomena:  the  evolution  of  the  EM 
field  over  time  as  it  approaches  a  steady  state,  and  the  steady  state 
solution.  As  described  in  Section  2.2  above,  our  system  computes 
EM  field  evolution  in  interactive  time  (several  updates  per  second), 
and  we  found  the  computational  steering  controls  to  be  an  effective 
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way  of  exploring  parameter  space  for  field  evolution.  However,  the 
steady  state  computation  is  on  the  order  of  minutes,  which  is  too 
slow  for  interactive  parameter  exploration.  This  motivated  the  de¬ 
velopment  of  the  objective  function  seeker  module  (see  Figure  1). 
This  module  allows  the  scientist  to  define  an  objective  function, 
which  is  a  formal  way  of  describing  the  particular  simulation  result 
which  is  desired.  An  example  of  an  objective  function  from  our 
application  area  is:  “ maximize  the  value  of  the  EM  field  intensity, 
for  all  points  inside  the  cylinder,  over  one  source  frequency  cycle ”, 
The  module  automatically  adjusts  steering  parameters  and  iterates, 
converging  on  the  objective  function.  This  module  uses  various  op¬ 
timization  techniques,  which  are  not  covered  in  this  paper,  but  are 
described  in  Smith  et  al.  [10]. 

3  Results  and  Lessons  Learned 

3.1  Visualizations 

The  microwave  simulations  produce  a  pair  of  vector  fields,  an  elec¬ 
trical  field  and  a  magnetic  field.  Visualizing  vector  fields  is  one  of 
the  most  challenging  problems  in  scientific  visualization,  because 
—  unlike  a  simple  scalar  field  of  amplitude  values  —  there  is  no 
direct  method  for  visualizing  a  vector  field.  Furthermore,  both  the 
field  domain  and  range  are  multi-dimensional,  and  so  understand¬ 
ing  the  underlying  topological  and  geometrical  structures  is  a  very 
challenging  task. 

Visualization  Techniques:  We  tried  a  number  of  techniques  to 
visualize  the  two  vector  fields.  In  particular,  we  tried  hedgehogs, 
streamlines,  particle  systems,  slice  planes,  and  volume  rendering. 
Each  of  these  methods  has  various  pros  and  cons. 

The  hedgehog  method  uses  arrow  icons  to  represent  the  3D  vec¬ 
tor  field  at  each  point  in  space.  The  length,  direction,  and  color  of 
the  arrow  can  visualize  all  the  parameters  of  a  3D  vector  field.  We 
found  that,  although  this  visualization  method  is  frequently  good 
for  understanding  global  structures,  it  does  not  easily  depict  local 
structures.  In  particular,  when  we  applied  hedgehogs  to  a  dense 
and  locally  varying  field  structure,  the  result  was  not  very  infor¬ 
mative  and  even  confusing,  due  to  the  clustering  and  overlapping 
of  vectors  in  a  small  region.  Streamlines  treat  the  3D  vector  as  a 
velocity  vector,  and  the  resulting  displacement  is  visualized  by  us¬ 
ing  continuous  lines  connecting  vectors  of  constant  velocity.  Since 
the  simulated  EM  fields  are  frequently  discontinuous  and  periodic, 
disjointed  and  repeating  patterns  of  the  streamlines  create  a  clus¬ 
tering  problem  similar  to  that  of  hedgehog  techniques.  A  particle 
system  is  a  variant  of  streamlines,  which  visually  animates  small 
particles  that  flow  along  the  path  of  the  displacement  described  by 
the  streamlines.  This  technique  provides  a  better  visualization  by 
avoiding  the  crowding  effect  found  in  both  hedgehogs  and  stream¬ 
lines,  but  we  found  that  it  was  not  suitable  for  visualizing  our  time- 
varying  vector  fields.  We  also  experimented  with  isosurfaces,  but 
we  found  them  unable  to  deal  well  with  the  complex  multimodal 
field  structures. 

Among  the  most  effective  were  simple  color-coded  slice  planes 
(see  the  lower  right-hand  window  in  Color  Plate  1),  which  did  a 
good  job  of  visualizing  field  intensity.  Since  the  more  interesting 
aspects  of  our  simulation  occur  in  specific  regions  —  for  exam¬ 
ple  along  the  center  axis  and  around  the  wall  slots  of  the  standard 
cylinder  —  we  can  place  a  slice  plane  through  these  regions  of  the 
visualization  space.  With  slice  planes  we  can  visualize  both  local 
and  global  field  features,  and  we  can  avoid  clustering  effects  by 
using  only  one  plane  for  visualization  through  a  region. 

Finally,  we  volume  rendered  the  field  energy  density  through¬ 
out  the  computational  space.  We  found  transparency  mapping  to 
be  the  key  parameter  for  creating  understandable  visualizations  — 
when  set  correctly,  internal  structures  were  easily  observable  inside 


the  field  volume.  Most  of  the  field  features,  both  global  and  local, 
can  be  displayed  by  volume  rendering.  While  there  are  many  ad¬ 
vantages  to  volume  rendering,  we  found  that  details  of  the  internal 
geometric  structure  of  the  field  were  difficult  to  examine  if  parts  of 
the  outer  layers  of  the  field  were  rendered  opaque  by  the  mapping 
function  that  we  used.  Color  Plate  4  shows  a  volume  rendering  of 
the  field  energy  density,  with  the  standard  cylinder  rendered  as  a 
wireframe  to  aid  in  visualization.  Note  the  two  hot  spots  which  can 
be  seen  in  the  field  on  each  side  of  one  of  the  slits  in  the  cylinder 
wall. 

Visualization  Issues:  As  mentioned  above,  our  simulation  actu¬ 
ally  produces  two  3D  vector  fields  —  an  electric  field,  and  a  mag¬ 
netic  field.  If  we  wish  to  visualize  both  fields  concurrently,  we  have 
to  utilize  two  different  techniques.  We  experimented  with  combina¬ 
tions  of  the  techniques  described  above.  In  particular,  we  had  good 
results  with  a  combination  of  hedgehogs  and  slicing  planes.  We 
found  that,  although  volume  rendering  is  good  for  visualizing  field 
density,  it  does  not  blend  well  with  other  visualization  techniques, 
because  it  fills  the  whole  computational  space  and  obscures  addi¬ 
tional  geometric  objects.  A  slice  plane,  on  the  other  hand,  occupies 
a  significantly  smaller  space,  and  thus  allows  other  geometric  ob¬ 
jects  to  be  seen. 

In  addition  to  the  3D  field  structure  at  a  given  time  step,  we  are 
also  interested  in  the  behavior  of  the  field  over  time.  Temporal  in¬ 
formation  comes  from  the  transient  behavior  of  the  field  structure 
as  it  evolves  from  time  step  0  to  a  steady  state.  At  steady  state  the 
EM  field  exhibits  a  stable  standing  wave  pattern,  which  is  easily 
seen  by  observing  the  field  at  every  nth  time  step,  where  n  is  the 
number  of  time  steps  per  cycle  of  the  EM  source;  at  steady  state  the 
pattern  is  stable  from  observation  to  observation.  Another  approach 
is  to  compute  and  display  the  average  magnitude  of  each  point  in 
the  field  every  n  time  steps.  This  presents  a  direct  view  of  the  time 
average  power  pattern  as  well  as  an  easy  way  to  assess  steady  state 
conditions. 

Virtual  Reality  Environments:  As  described  in  Section  2.3 
above,  in  addition  to  viewing  our  results  via  the  AVS  geometry 
viewer  on  our  computer  monitor,  we  used  our  VR  Viewer  mod¬ 
ule  to  view  the  results  in  our  CAVE  facility  (Color  Plate  3).  We 
found  several  advantages  here.  First,  since  the  CAVE  screens  are 
much  larger  than  normal  computer  monitors,  we  were  able  to  see 
images  that  completely  surround  us  with  a  very  large  field  of  view. 
In  addition,  we  could  render  the  visualization  in  stereo,  which  made 
3D  effects  more  prominent.  We  could  interact  with  a  hand-held  3D 
wand,  which  allowed  us  to  probe  and  manipulate  the  field  images 
directly.  We  believe  the  CAVE  environment  helped  us  to  better  un¬ 
derstand  our  data. 

We  also  found  advantages  to  driving  the  CAVE  display  from  our 
VR  Viewer.  The  main  advantage  is  that  the  VR  rendered  image  is 
part  of  an  AVS  visualization  network.  Thus  we  can  interactively 
perform  all  the  visualization  functions  provided  by  AVS,  as  well  as 
our  other  custom  modules.  The  same  advantage  extends  to  end-user 
scientists  as  well:  they  can  visualize  their  scientific  information  in  a 
virtual  environment  (VE)  with  the  same  visualization  toolkit  (AVS) 
that  they  are  already  using.  In  addition,  VE  environments  such  as 
CAVEs  are  still  relatively  uncommon,  and  may  not  be  located  near 
a  scientist’s  office.  The  VR  Viewer  allows  scientists  to  develop  their 
visualization  networks  at  their  own  desktops,  and  only  visit  the  VE 
facility  when  they  have  data  to  visualize. 

Automatic  Viewing  Parameter  Tuning:  We  have  experi¬ 
mented  with  automatically  tuning  viewing  parameters  for  complex 
visualization  applications.  In  particular,  we  have  developed  a  mod¬ 
ule  which  tunes  hedgehog  displays  of  vector  fields  [4],  If  the  scale 
of  the  vector  arrows  is  too  small,  or  the  density  of  arrows  is  too  low, 
the  vector  information  may  not  appear  in  the  display.  Conversely,  if 
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the  scale  is  too  large,  or  the  density  of  arrows  is  too  high,  then  the 
display  will  be  saturated  with  arrows  and  all  other  information  will 
be  obliterated.  Without  the  automated  parameter  selection,  the  user 
must  test  various  parameter  values  by  hand. 

VolumePro  Visualization:  In  a  recent  development,  we  ob¬ 
tained  a  VolumePro  real-time  volume  rendering  board  [8],  which 
fits  in  a  high-end  PC.  We  have  begun  experimenting  with  visualiz¬ 
ing  our  data  using  this  system;  Color  Plate  5  shows  a  preliminary 
result.  The  advantages  of  the  VolumePro  board  are  that  it  provides 
high-quality  volume  rendering  at  very  interactive  rates,  allows  in¬ 
teractive  rendering  not  only  of  still  volumes,  but  also  of  animated 
volumes,  and  of  course  that  it  costs  an  order  of  magnitude  less  than 
an  equivalent  high-end  graphics  computer.  We  have  several  future 
plans  for  the  VolumePro  board  in  this  project.  First,  we  are  devel¬ 
oping  software  that  will  provide  real-time,  mouse-based  updating 
of  color  maps,  while  the  volume  dataset  is  being  animated.  As  find¬ 
ing  the  color  map  that  best  displays  the  relevant  information  can 
be  quite  difficult  and  time-consuming,  this  should  speed  up  this  as¬ 
pect  of  using  a  visualization  system.  Second,  we  hope  to  integrate 
the  VolumePro-enabled  PC  into  our  overall  computation  steering 
system.  This  will  involve  writing  another  custom  AVS  module, 
which  will  drive  a  volume  rendering  application  on  the  VolumePro- 
enabled  PC. 

3.2  Inverse  Steering 

We  implemented  the  inverse  steering  module  discussed  in  Sec¬ 
tion  2.4  above,  and  used  it  to  explore  a  number  of  objective 
functions.  However,  because  our  supercomputer  required  about  a 
minute  to  achieve  a  steady  state  condition  for  each  steering  param¬ 
eter  tuple,  even  an  optimal  search  algorithm  required  several  hours 
to  find  a  solution  over  the  full  parameter  5-space.  In  addition,  we 
found  it  common  for  this  solution  to  exhibit  unforeseen  and  undesir¬ 
able  characteristics.  Although  we  could  refine  the  objective  func¬ 
tion  to  avoid  solutions  with  these  characteristics,  at  several  hours 
per  iteration  we  found  ourselves  back  in  the  realm  of  batch  com¬ 
puting.  We  wanted  to  find  a  way  to  make  the  inverse  steering  more 
interactive. 

We  therefore  performed  the  following  experiment.  We  pre¬ 
computed  the  steady  state  field  result  for  22,638  parameter  points 
uniformly  sampled  from  the  4-space  of  source  frequency  (33  sam¬ 
ples),  source  direction  angle  (7  x  14  samples),  and  polarization 
angle  (7  samples).  We  left  out  probe  insertion  distance  to  make  the 
problem  tractable,  and  because  we  are  most  interested  in  the  no¬ 
probe  condition.  This  computation  took  about  25  hours  of  super¬ 
computer  time  and  resulted  in  a  6  gigabyte  database.  We  developed 
a  simple  interface  to  search  the  database,  and  return  the  field  result 
for  any  particular  4-tuple  of  steering  parameters.  This  gives  us  a 
simulated  HASP  computation,  which  takes  only  about  10  millisec¬ 
onds  to  return  a  steady  state  field  result.  Using  this  database,  we 
can  run  the  inverse  steering  module  interactively.  This  gives  us  all 
the  usual  benefits  of  interactivity  as  compared  to  batch  processing, 
although  of  course  we  are  limited  to  the  particular  quantization  of 
the  parameter  values  and  field  results  which  we  used  to  create  the 
database. 

An  interesting  result  from  this  pre-computed  database  is  that  we 
can  compute  any  objective  function  over  the  entire  parameter  space, 
and  visualize  the  results  using  our  standard,  forward  steering  inter¬ 
face.  Color  Plate  2  is  an  example  of  this.  Here  we  are  visualizing 
the  peak  EM  power  intensity  (the  same  objective  function  specified 
in  Section  2.4  above),  as  a  function  of  three  of  the  five  steering  pa¬ 
rameters.  Notice  in  Color  Plate  2a  that,  in  addition  to  local  “hot 
spots”  representing  higher  values  in  the  parameter  space,  there  is  a 
global  bimodal  structure  with  a  “cloud”  of  higher  peak  power  val¬ 
ues  occurring  in  the  lower  frequency  region  (around  1  GHz),  and 
a  second  smaller,  less  intense  cloud  of  peak  power  values  in  the 


higher  frequency  region  (around  2  GHz).  This  bimodal  structure  is 
consistent  with  the  analytic  observation  that  the  standing  wave  pat¬ 
tern  in  a  closed  cylinder  would  tend  to  peak  at  wavelengths  which 
are  submultiples  of  each  other. 

4  Conclusions 

In  this  paper,  we  have  described  a  computational  steering  system 
architecture,  which  we  developed  to  support  the  execution  and  visu¬ 
alization  of  microwave  interactions  with  missile  bodies.  We  found 
the  primary  advantages  of  our  system  were  that  1)  it  is  interactive, 
which  allows  visualization  parameters  to  be  more  thoroughly  ex¬ 
plored,  2)  its  interface  is  implemented  within  the  standard  visual¬ 
ization  package  AVS,  3)  it  allows  scientists  to  view  and  manipulate 
their  visual  output  in  a  large-format  CAVE  environment,  and  4)  it 
supports  inverse  steering,  which  enhances  the  ability  of  scientists  to 
find  simulation  parameters  which  produce  desired  results. 
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Color  Plate  2:  (a)  Volume  rendering  of  the  standard 
cylinder  peak  power  intensity  as  a  function  of  three 
steering  parameters:  EM  source  frequency,  EM  source 
location  (0,  4»).  (b)  Slices  through  the  two  resonant 

regions  at  1.0  and  2.0  GHz. 


Color  Plate  4:  Volume  rendering  of  the  field  together 
with  a  wireframe  rendering  of  the  standard  cylinder. 


Color  Plate  5:  The  dataset  rendered  with  the  VolumePro 
board. 


